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ABSTRACT 



The spectrum from the black hole X-ray transient GRO J1655-40. obtained using the 
Chandra High Energy Transmission Grating (HETG) in 2005 is notable as a laboratory 
for the study of warm absorbers, and for the presence of many lines from odd-Z elements 
between Na and Co (and Ti and Cr) not previously observed in X-rays. We present 
synthetic spectral models which can be used to constrain these element abundances 
and other parameters describing the outflow from the warm absorber in this object. We 
present results of fitting to the spectrum using various tools and techniques, including 
automated line fitting, phenomenological models, and photoionization modeling. We 
show that the behavior of the curves of growth of lines from H-like and Li-like ions 
indicate that the lines are either saturated or affected by filling-in from scattered or 
a partially covered continuum source. We confirm the conclusion of previous work by 



Miller et al. (2006) and Miller et al. (20081 which shows that the ionization conditions 



are not consistent with wind driving due to thermal expansion. The spectrum provides 
the opportunity to measure abundances for several elements not typically observable 
in the X-ray band. These show a pattern of enhancement for iron peak elements, and 
solar or sub-solar values for elements lighter than calcium. Models show that this is 
consistent with enrichment by a core-collapse supernova. We discuss the implications 
of these values for the evolutionary history of this system. 
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Introduction. 



X-ray spectra reveal that warm absorbers (absorption by partially ionized gas) are a common 
feature in compact objects. Although warm absorbers were first detected in the spectra of active 
galaxies using the Einstein satellite ( Halpern|1984 |, Chandra and XMM — Newton grating obser- 
vations have shown that these occur in X-ray binaries as well, and information about the dynamics 
and other properties of this gas can have important implications for our understanding of accretion 
in these systems. A notable example is the 900 ks observation using the Chandra High Energy 



Transmission Grating (HETG) of the Seyfert galaxy NGC 3783 ( |Kaspi et ar]|2002[ ), but this has 
been surpassed in signal-to-noise by the 2005 spectrum of the galactic black hole transient GRO 
J1655-40 (Miller et al. 20061. The statistical quality of this spectrum makes it one of the best 



available for testing and refinement of models for warm absorber flows. 

GRO J 1655-40 was discovered using BATSE onboard the Compton Gamma Ray Observatory 



(Harmon et al. 19951. Radio observations show apparently superluminal jets (Hjellming & Rupen 



19951. Optical observations obtained during quiescence show the companion star is a F3-5 giant 



or sub-giant in a 2.62 day orbit around a 5-8 M© compact object (Orosz Sz Bailyn 1997 1 . Deep 



absorption dips have been observed ( Balucihska-Church|200l" | suggesting inclination ~ 70° ( [van der 



Hooft et al. 1998). GRO J1655-40 shows the highest- frequency quasi-periodic oscillations (QPOs) 
seen in a black hole ( Strohmayer||200T I . Narrow X-ray absorption lines from highly ionized Fe were 



detected by Yamaoka et al. (2001 1 and Ueda et al. (19981. Similar features were detected from all 



the bright dipping low mass X-ray binaries (LMXBs) observed with XMM- Newton (Diaz Trigo et 
al] |2006p , and from the LMXB GX 13-^1 ( |Ueda et al!]|2001[ ). This suggests that ionized absorbers 
are a common feature of LMXBs, although they may not be detected in objects viewed at lower 
inclination. 



The properties of the GRO J1655-40 warm absorber have been explored by Miller et al. (2006 1 



and Miller et al. (20081. The spectrum resembles warm absorbers observed from other compact 



objects in that the lines are blueshifted, and that the inferred Doppler blueshift velocities are in 
the range 300 - 1600 km/s. The lines are identified primarily as H- and He-like species of elements 
with nuclear charge 8 < Z < 28, and there is no clear evidence for absorption by low ionization 
material such as the iron M-shell UTA (Behar et al. 2003). The unprecedented signal-to-noise 



may account for the presence of lines from many trace elements previously not detected in X-rays, 
essentially all elements between Na and Co. The detection of the 11.92 A line, arising from the 
2^3/2 metastable level of Fe XXII, implies a gas density > 10^^ cm~^. The spectrum is richer in 
features than another spectrum taken earlier in the same outburst, which showed only absorption 
by Fe XXVI La. The reason for this richness may be due to the higher column density of absorber 
and to the softer continuum during this particular observation, although this has not been tested 
quantitatively. 



The observed line strengths can be used to infer the ionization balance, i.e. the ratio of 
abundances of H-like, He-like and lower ion stages for various elements. Models for the ionization 
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balance in the wind then yield the ionization parameter: the ratio of the ionizing flux to the gas 
density, and the density is constrained by the Fe XXII line detection. This, together with the 



observed luminosity and the observed outflow speed, led Miller et al. (20061 to conclude that the 



radius where the outflow originates is too small to allow a wind driven by thermal pressure. That is, 
the likely ion thermal speeds in the gas are less than the escape velocity at the inferred radius. On 
the other hand, the mass flux can be inferred from the line strengths, and the estimated rate exceeds 
what is expected from outflows driven by radiation pressure. This implies that the outflow must 
be driven by a different mechanism, such as magnetic forces, but this result has been controversial 



(Netzer 20061. 



The mass flux in the wind is important for our understanding of the mass and energy budget 
in accreting compact objects and it is clear that accurate models for the ionization balance and 
synthetic spectrum are needed in order to reliably determine the properties of the wind. It is the 



goal of this paper to study this, and several issues which were not considered by Miller et al. (2006 



2008): (1) What element abundances are required to account for the lines from the many iron 



peak elements observed in the spectrum, and what might this tell us about the origin of the gas 
in the wind and accretion flow? This is the only spectrum obtained from a warm absorber which 
clearly detects lines from odd-Z elements with Z > 10 and from iron peak elements other than Fe 
and Ni. The abundances are of interest since they may contain clues to the evolutionary origin of 



the system; Israelian et al. (1999) find evidence for enhanced 0/H, Mg/H, Si/H, S/H, and relative 



to solar, but not Fe/H, and suggest that this may be due to enrichment by the supernova which 
produced the compact object. (2) What is the possible influence of radiative transfer effects (line 
emission, partial covering, or flnite energy resolution) on the inferred wind properties? Analyses 
of X-ray warm absorber flows so far have not attempted to account in detail for these processes 



(although hints to their importance in AGN come from the UV; Gabel et al. (2005)). They could 
systematically skew the derived column densities and ionization conditions. It is straightforward, 
although time consuming, to construct models which will test these effects in various scenarios. (3) 
Are there correlations between line shape or centroid and the ionization conditions where that line 
is expected to dominate? A flow with a predominantly ordered velocity fleld and a central radiation 
source is likely to show a gradient in ionization balance with position, and therefore with velocity. 
This should be manifest as correlations between line width or offset with the ionization degree of 



its parent ion. This effect is not found in Seyfert galaxy warm absorbers (Kaspi et al. 2002), but 



examination of the GRO J 1655-40 spectrum shows differences among the line profile shapes. After 
careful attention to item (2) above, we can test this quantitatively. 

In the remainder of this paper we present our model fits and interpretation of the Chandra 
HETG spectrum of GRO J1655-40, attempting to address the above questions. This includes 
various fitting techniques, testing of radiative transfer effects, and discussion of the constraints on 
element abundances. Finally, we discuss the implications of these results for the dynamics of the 
outfiow, and possible evolutionary history of the source. 
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2. Fitting: Notch Models 



The fits in this paper were performed using the same extraction of the HETG data as was 



used by Miller et al. (20081. Chandra observed GRO J1655-40 for a total of 44.6 ks starting 



at 12:41:44 UT on 2005 April 1. Data was taken from the ACIS-S array dispersed by the High 
Energy Transmission Grating Spectrometer (HETGS). Continuous clocking mode was used in order 



to prevent photon pile-up. As described by Miller et al. (2008), a gray filter was used in order to 



reduce the zero order counting rate. Data was processed using the CIAO reduction package, version 
3.2.2. The event file was filtered to accept only standard event grades, good-time intervals, and 
to reject bad pixels. Streaking was removed using the "destreak" tool, and spectra were extracted 
using "tg jesolve_events" and "tg_extract" . Arfs were produced using the "fullgarf" tool along with 
canned rmfs. First order HEG spectra and arfs and first-order MEG spectra and arfs were added 
using the "add_grating^pectra" tool. In this paper we do not fit to the RXTE data obtained 
simultaneously with the Chandra HETG observation, but we do employ the continuum shape 



derived from the Miller et al. (20081 fits which include the RXTE data. 



Our procedure for analyzing the spectrum of GRO J1655-40 consists of three separate parts. 
First, we fit the spectrum using xspec together with simple analytic models describing the contin- 
uum and the lines. This consists of a continuum shape which is the same as that used by [Miller et 



al. (2008): a power law plus a disk blackbody and cold absorption. Based on the fits by Miller et al 



(20081, we use a power law index of 3.54 and a disk blackbody temperature with temperature 1.34 



keV. We allow the normalizations of the components to vary, and find a best fit normalization of 
515.7ibl.5 for the disk blackbody and <0.15 for the power law. The flux is 2.02 xlO~® ergs/cm^/s 
2-10 keV. We ignore photons with wavelength greater than 15 A, since there are few counts in this 
range and its inclusion has negligible effect on the results at shorter wavelength. We refer to this 
as model 1. The fitting parameters and for this and the other fits discussed in this paper are 
presented in Table [T] The various physical parameters for all the models we test are listed in the 
first column. Model 1, the best fit to the continuum only, gives x^/'^ = 120618/8189. 
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Table 1. Model results 



fitting quantity 


units 


Model 1 


Model 2 


Model 3 


Model 4 


Model 5 


Model 6 


log(«l) 


erg cm 


NA 


NA 


NA 


4.I0.I 


s.slo.i 


4.0I0.I 


log(Ni) 


_9 

cm 


NA 


NA 


NA 


23.8l0.02 


22.64l0.02 


24.0l0.02 




km 


NA 


NA 


NA 


50* 


200* 


200* 




km 


NA 


NA 


NA 


-375 


-375 


-375 


iog(6) 


erg cm 


NA 


NA 


NA 


NA 


4.6I0.I 


NA 


log(N2) 


cm~^ 


NA 


NA 


NA 


NA 


23.90 


NA 


EW/e26 


keV 


NA 


NA 


NA 


0.03 


NA 


0.03 


Vo//,/e26 


km 


NA 


NA 


NA 


1451. 


NA 


1451. 


Vtur(),/e26 


km 


NA 


NA 


NA 


n 03+°°^ 


NA 


n 03+°'^^ 


NH 


cm~^ 


21.49l0.01 


21.49l0.01 


NA 


21.49l0.01 


21.48l0.01 


21.48l0.01 


7 




3.54 


3.54 


3.54 


3.54* 


3.54* 


3.54* 


pi norm 




<0.01 


<0.01 


<0.15 


<0.01 


<0.01 


<0.01 


diskbb norm 




513ll 


538* 


534+ _1 


548l?-J 




533li-^ 




keV 


1.35 


1.35* 


1.35* 


1.35* 


1.35* 


1.35* 


fscatt 



















0.37 






120618 


18474 


20671 


33591 


34283 


24561 


dof 




8189 


8108 


8118 


8187 


8189 


8184 



* These parameters were fixed during the fitting. 
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We then add the effects of absorption hnes by using negative Gaussian models for the hnes, as 



was done by Miher et ah (20081. The hst of hnes and their properties is the primary topic of the 



rest of this paper. InitiaUy we use the same hnes as given in Table 1 of Miller et al. (20081. There 
are 71 lines with well-determined wavelengths and identifications in this list. We allow the values 
of the wavelengths, widths, and line normalizations to be determined by the XSPEC minimization. 



This results in detections of essentially all the lines with parameters consistent with those of Miller 



et al. (2008). In addition, we propose identifications for some of the lines which were not identified 
by those authors. We point out that this procedure differs from those authors in that we fit to a 
single analytic global continuum, while they fit to a piece-wise powerlaw. The procedure used here 
is chosen for comparison with the fits to photoionization models in the next section. The best-fit 
to the continuum plus Gaussians gives 'x^ jv =18474/8108. We refer to this as model 2 (cf. Table 
|l). This fit is marginally acceptable based on standard arguments derived from statistics, and 
it is the best of the models presented in this paper. This serves to illustrate the level of systematic 
errors, or errors in our continuum, which provide an effective limit to our ability to fit the spectrum. 
Model 2, along with other models discussed below, is plotted in figures [T]- 114[ In this figure the 
vertical axis is the ratio of the model and data to the continuum-only model, model 1. Successive 
models are offset by unity from each other. These figures show qualitatively the agreement between 
the model and the data which we achieve. 

The second approach to line fitting is to construct an automated line fitting program. This 
takes the same continuum employed in the continuum-only fit, and then experiments with random 
placements of lines throughout the 1-15 A range. These experiments begin with an initial wave- 
length which is chosen randomly within this range (but excluding a region within 2 Doppler widths 
of previously found lines), and then the line wavelength, width, and optical depth are varied in an 
attempt to find a best fit (the wavelength is restricted to a region near the initial wavelength in 
this procedure). The fit is considered valid if the improves by 3 with the inclusion of the line. 
Line widths are limited to be less than 8 Doppler widths when compared with a turbulent velocity 
of 100 km/s. The lines are assumed to be Doppler broadened only, and the absorption is a true 
Gaussian absorption profile. This is in contrast to the xspec Gaussian line model, which treats 
absorption as a negative emission. As we will show, many of our line fits result in large optical 
depths, and in this case a Voigt profile fit would be preferable. The pure Doppler profile likely 
results in an over-estimate of the line optical depth, since it cannot produce as strong absorption 
line wings as would a Voigt profile. However, the Voigt profile damping parameter value depends 
on the line identification, and cannot be conveniently used as a fitting parameter. In our synthetic 
spectral modeling, in the following section, we fit to Voigt profiles using accurate atomic rates for 
the damping parameters. For convenience we refer to this as the notch model, although it assumes 
Gaussian absorption lines rather than true notches. 

This procedure yields a total of 292 lines after a total of 20000 attempts at placing random 
lines. We neglect lines with equivalent widths less than 3.4 xlO~^ eV. We do not consider this to 
be necessarily an exhaustive list of lines in the GRO J 1655-40 spectrum, but likely contains the 
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strongest or least ambiguous features, and it is objective. This procedure detects all the features 
in the Miller et al. (20081 table, plus many others which are weaker or blended. Some of these 
are undoubtedly artifacts of the shortcomings of our continuum model, and others may coincide 
with regions where bound-free continuum opacity is important. These are discussed in turn in the 
following. This fit yields ^jv = 20671/8118 using the continuum from the XSPEC Gaussian fit. We 
refer to this as model 3 (cf. Table [l|. The slightly worse value for this model when compared 
with model 2 is likely due to the limitations of our automated fitting procedure, particularly when 
two strong lines are close together or partially overlapping. In addition to deriving wavelengths, line 
center optical depths, and widths (a), we also derive errors on these quantities based on the Ax^=3 
criterion of Cash ( 1979 1, and we calculate the line equivalent widths by integrating numerically over 
the best-fit model profile. 



2.1. Line Identifications 

The list of lines we detect is given in Table |2] This includes the wavelength, width (a) 
and equivalent width derived from the automated fitting leading to model 3. We also provide 



identifications for the lines. This is done by searching the linelist in the XSTAR (Kallman and 



Bautista||2001 Bautista and Kallman|[2001 ) database and choosing the line which has the greatest 



ratio of optical depth to Doppler shift within a Doppler shift of <ll500 km s~^. The optical 
depth used in this determination is calculated using the WARMABS analytic model, which is an 
implementation of xstar use as an analytic model within the xspec X-ray spectral fitting package, 
for the conditions in model 4 described in the following section. We have updated both xstar and 
WARMABS to include all the previously neglected elements with Z < 30; this is described in the 
Appendix. Note that the criterion for line identification is used to choose among the xstar lines 
which fall within Doppler shifts of ^'^1500 km s~^, but does not prevent a line from being included if 
there is no xstar line within that interval. If there is an ID, then the parent ion, xstar wavelength, 
and upper and lower level designations are also given in Table |2] The identification, together with 
the optical depth derived from the model 3 fits allows the equivalent hydrogen column density of 
the absorbing ion to be derived. These are discussed in more detail in the following subsection, and 
are given in the table. The elemental abundances used in the calculation of equivalent hydrogen 
abundances are those of ( Grevesse et al.|1996 Allen|1973 1. Figures [I] - [l4| show the count spectrum 



observed by the HETG (relative to the continuum only model 1) together with the lines and 
identifications from Table |2j Also shown in these figures is our fit to model 2 and to models 4, 5 
and 6 which will be discussed in more detail below. 
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Fig. 1. — Spectrum AA 1 - 2 A. Spectrum is shown as ratio relative to pure power law model (model 
1). Various models 2,4,5,6 are labeled. The vertical axis is the ratio of the model and data to the 
continuum-only model, model 1. Successive models are offset by unity from each other. 
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Fig. 2. — Spectrum AA 2 - 3 A. Spectrum is shown as ratio relative to pure power law model (model 
1). Various models 2,4,5,6 are labeled. The vertical axis is the ratio of the model and data to the 
continuum-only model, model 1. Successive models are offset by unity from each other. 
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Fig. 3. — Spectrum AA 3 - 4 A. Spectrum is shown as ratio relative to pure power law model (model 
1). Various models 2,4,5,6 are labeled. The vertical axis is the ratio of the model and data to the 
continuum-only model, model 1. Successive models are offset by unity from each other. 
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Fig. 4. — Spectrum AA 4 - 5 A. Spectrum is shown as ratio relative to pure power law model (model 
1). Various models 2,4,5,6 are labeled. The vertical axis is the ratio of the model and data to the 
continuum-only model, model 1. Successive models are offset by unity from each other. 
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Fig. 5. — Spectrum AA 5 - 6 A. Spectrum is shown as ratio relative to pure power law model (model 
1). Various models 2,4,5,6 are labeled. The vertical axis is the ratio of the model and data to the 
continuum-only model, model 1. Successive models are offset by unity from each other. 
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Fig. 6. — Spectrum AA 6 - 7 A. Spectrum is shown as ratio relative to pure power law model (model 
1). Various models 2,4,5,6 are labeled. The vertical axis is the ratio of the model and data to the 
continuum-only model, model 1. Successive models are offset by unity from each other. 
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Fig. 7. — Spectrum AA 7 - 8 A. Spectrum is shown as ratio relative to pure power law model (model 
1). Various models 2,4,5,6 are labeled. The vertical axis is the ratio of the model and data to the 
continuum-only model, model 1. Successive models are offset by unity from each other. 
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Fig. 8. — Spectrum AA 8 - 9 A. Spectrum is shown as ratio relative to pure power law model (model 
1). Various models 2,4,5,6 are labeled. The vertical axis is the ratio of the model and data to the 
continuum-only model, model 1. Successive models are offset by unity from each other. 
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Fig. 9. — Spectrum AA 9 - 10 A. Spectrum is shown as ratio relative to pure power law model 
(model 1). Various models 2,4,5,6 are labeled. The vertical axis is the ratio of the model and data 
to the continuum-only model, model 1. Successive models are offset by unity from each other. 
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Fig. 10. — Spectrum AA 10 - 11 A. Spectrum is shown as ratio relative to pure power law model 
(model 1). Various models 2,4,5,6 are labeled. The vertical axis is the ratio of the model and data 
to the continuum-only model, model 1. Successive models are offset by unity from each other. 
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Fig. 11. — Spectrum AA 11 - 12 A. Spectrum is shown as ratio relative to pure power law model 
(model 1). Various models 2,4,5,6 are labeled. The vertical axis is the ratio of the model and data 
to the continuum-only model, model 1. Successive models are offset by unity from each other. 
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Fig. 12. — Spectrum AA 12 - 13 A. Spectrum is shown as ratio relative to pure power law model 
(model 1). Various models 2,4,5,6 are labeled. The vertical axis is the ratio of the model and data 
to the continuum-only model, model 1. Successive models are offset by unity from each other. 
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Fig. 13. — Spectrum AA 13 - 14 A. Spectrum is shown as ratio relative to pure power law model 
(model 1). Various models 2,4,5,6 are labeled. The vertical axis is the ratio of the model and data 
to the continuum-only model, model 1. Successive models are offset by unity from each other. 
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Fig. 14. — Spectrum AA 14 - 15 A. Spectrum is shown as ratio relative to pure power law model 
(model 1). Various models 2,4,5,6 are labeled. The vertical axis is the ratio of the model and data 
to the continuum-only model, model 1. Successive models are offset by unity from each other. 
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Table 2. Line Properties 





index 




A ' 


ion 






ew 


width* 


1.544 


17264 


1 


.542 


Mn XXV 


-357, 


,5 


0, 


,0006 


459.3 


1.567 


128239 


1 


.573 


Fe XXV 


1187, 


,0 


0, 


,0029 


772.1 


1.582 


132917 


1 


.589 


Ni XXVII 


1232, 


,6 





,0084 


506.7 


1.601 








no ID 






0, 


.0022 


707.2 


1.613 


17063 


1 


.615 


Cr XXIV 


297, 


,6 





.0004 


234.9 


1.626 


17263 


1 


.626 


Mn XXV 


77, 


,5 


0, 


,0010 


412.7 


1.646 


128543 


1 


.649 


Co XXVII 


484, 


,8 


0, 


,0006 


407.0 


1.675 


17061 


1 


.674 


Cr XXIV 


-107, 


,5 


0, 


,0004 


204.0 


1.693 








no ID 






0, 


.0007 


458.1 


1.709 


128440 


1 


.712 


Co XXVI 


544, 


,2 





.0015 


539.4 


1.726 








no ID 









,0012 


701.0 


1.742 








no ID 






0, 


,0003 


225.4 


1.772 


128358 


1 


.780 


Fe XXVI 


1325, 


,6 


0, 


,0053 


453.6 


1.838 








no ID 






0, 


.0007 


195.2 


1.851 


128293 


1 


.850 


Fe XXV 


-81, 


,0 





.0157 


503.9 


1.865 


127741 


1 


.864 


Fc XXIV 


-209, 


1 


0, 


,0006 


179.5 


1.875 


127760 


1 


.873 


Fc XXIV 


-320, 


,0 





,0003 


203.7 


1.923 


17262 


1 


.926 


Mn XXV 


480, 


,5 





,0007 


401.6 


1.951 


16659 


1, 


.948 


Ti XXII 


-412, 


,1 


0, 


.0002 


299.7 


1.963 








no ID 









.0003 


212.7 


1.993 


16658 


1 


.995 


Ti XXII 


313, 


1 





,0003 


207.6 


2.004 


17159 


2 


.006 


Mn XXIV 


329, 


,4 


0, 


,0018 


566.0 


2.020 








no ID 






0, 


,0007 


493.5 


2.036 








no ID 






0, 


.0003 


454.9 


2.062 








no ID 









.0003 


442.3 


2.088 


17059 


2 


.092 


Cr XXIV 


502, 


,9 


0, 


.0014 


416.2 


2.114 








no ID 






0, 


.0002 


161.4 


2.153 








no ID 






0, 


.0005 


156.3 


2.179 


16956 


2, 


.182 


Cr XXIII 


426, 


,8 


0, 


.0020 


526.2 


2.193 


17016 


2 


.193 


Cr XXIII 


-68, 


,4 





.0004 


181.7 


2.205 








no ID 






0, 


.0002 


158.5 


2.312 








no ID 









.0003 


158.4 


2.361 








no ID 






0, 


,0003 


396.3 


2.417 








no ID 






0, 


.0004 


403.5 


2.491 


16656 


2 


.492 


Ti XXII 


150, 


,5 





.0007 


410.0 


2.502 


16183 


2 


.514 


Ca XIX 


1414, 


,8 


0, 


.0003 


151.4 


2.545 


16282 


2 


.549 


Ca XX 


515, 


,1 


0, 


.0011 


136.9 


2.555 








no ID 






0, 


,0002 


140.0 


2.701 


16143 


2, 


.705 


Ca XIX 


488, 


,7 


0, 


.0008 


400.9 


2.857 








no ID 









.0003 


117.3 


2.877 








no ID 






0, 


.0003 


124.1 


2.982 


13532 


2 


.987 


Ar XVIII 


548, 


,3 


0, 


.0007 


131.4 


3.016 


16262 


3 


.020 


Ca XX 


427, 


,7 


0, 


.0050 


415.2 


3.145 


13528 


3 


.151 


Ar XVIII 


534, 


,2 


0, 


.0016 


304.3 


3.172 


16146 


3 


.177 


Ca XIX 


491, 


,8 





.0024 


203.9 


3.355 


13463 


3 


.365 


Ar XVII 


938, 


,9 


0, 


.0009 


476.4 



lower level upper level ion column* 



Q 71+0. ,56 
^'•'-'--0.41 

23.06l|-?» 

1 10+0. 80 

9+0. .SI 
^-0.40 
fi+0.34 
-0.28 

ri + 0.79 

"^-0.63 
0+0.27 
-0.23 
, 24+0.40 

fi5+''-24 

U.00_o 22 

10.56l|?« 

9 QK+0.75 

0+9.99 

-8.98 
.,+0.96 



1.227 

0.66^ 
1.22^ 
0.58^ 



39.38^ 



-0.70 
+ 0.28 
0.24 
+ 0.13 
-0.12 
+0.14 
0.09 
+0.23 
-0.20 
.+0.21 
0.19 
0.15 
-0.14 
0.10 
0.09 
.+0.09 
0.09 
+0.08 
0.08 

0.84l°;13 

26+°-20 
'J-^°-0.17 

n 7Q+0.28 
'^•'•^-0.26 

0.87l«;ll 
0.52+0-20 



0.41] 
0.44^ 
0.14^ 
0.44^ 
0.36^ 
0.87^ 
0.32^ 
0.17^ 
O.M"' 



-0.18 
n+0.20 
^-0.14 

0.27+0-19 



0.22"' 



-0.17 
n+0.06 
-0.05 

0.15+0-06 



O.IO"! 



-0.05 
n+0.05 
^-0.05 

0.30+0-16 



0.22"' 



i.eo"! 



-0.14 
n+0.29 
-0.26 

20+0-1^ 
'^•^'^-0.13 

21+0-05 
"•^1-0.05 

0.18+0-13 



-0.11 
Q+0.16 
-0.15 

0.48+0-15 



0.23"! 



-0.13 
4+0.08 
^-0.07 

0.40+0-06 



1.24"* 



-0.05 

0+0.12 

-0.12 
11+003 



1.08"! 



0, 


,029 


ls.2S 


4p.2P° 


2.9+2-3 


X 


1023 


0, 


,152 


ls2.1S 


ls.3p.lP* 




X 


1021 


0, 


,684 


ls2.1S 


ls.2p.lP 




X 


1023 



0.008 


ls.2S 


6p.2P° 


1 4+0.9 
l-*-0.7 


X 


1024 


0.080 


ls.2S 


3p.2P° 


1 7+0.7 
1- ' -0.6 


X 


1023 


0.422 


ls.2S 


2p.2P° 


7 4+3.8 
'•^-3.1 


X 


1022 


0.029 


ls.2S 


4p.2P° 


7+2.4 
■^■'-1.9 


X 


1023 



0.693 



0.408 



0.029 
0.711 



0.421 



0.721 
0.152 



0.152 



ls2.1S 



ls.2S 



ls.2p.lP 8.2 



+2.6 
2.1 



X 10^ 



2p.2P° 5.912-1 X 10' 



0.775 


ls2.1S 


ls.2p.lP* 


1, 


1+0.3 
■^-0.3 


X 


1022 


0.149 


ls2.2s 


ls2s2p.2P0.5 


3, 


7+1.4 

■ ' -1.0 


X 


1021 


0.015 


ls2.2s 


ls2s2p.4P1.5 


5, 


7+3.9 
■'-3.4 


X 


1021 


0.420 


ls.2S 


2p.2P° 


10, 


0+2-9 
,U_2 7 


X 


1021 


0.014 


ls.2S 


5p.2P° 


4, 


1+4.0 
■1-2.6 


X 


1023 



ls.2S 
ls2.1S 



ls.2S 



ls2.1S 
ls2.1S 



4p.2P° 
ls.2p.lP 



2p.2P° 



ls.2p.lP 
ls.2p.3P 



5.0I2-? X 10^ 



1.1 



+0.2 
0.2 



1 4+0-2 



X 10-^ 



X 10^ 



8.2li-0 X 

9 ^i+0.9 V 



1021 
1022 



0.419 


ls.2S 


2p.2P° 


1 7+0.4 


X 


1022 


0.027 


ls2.1S 


ls.5p.lP* 


9 5+1-3 


X 


1022 


0.078 


ls.2S 


3p.2P° 


4 8+0-8 
+-°-0.7 


X 


1022 



ls2.1S 



ls.3p.lP* 



9 Q+0.7 



1021 



0.029 


ls.2S 


4p.2P° 


-I c+O.S 
l-''-0.4 


X 


1022 


0.411 


ls.2S 


2p.2P° 


5 7+0-4 
'^•'-0.3 


X 


1021 


0.078 


ls.2S 


3p.2P° 


4 3+0-6 


X 


1021 


0.770 


ls2.1S 


ls.2p.lP* 


9 c+O.S 
^■^-0.3 


X 


1021 


0.153 


ls2.1S 


ls.3p.lP* 


5.511:1 


X 


1020 
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index 



width<= 



o.oyu 


1 1 Q'Tfi 

iio/ O 


Q fiflfi 


o AVi 


A'7Ci '7 

4 f y. ( 


o noi o 
U.UUiz 


1 A'7 
iU / . 


.i 


A 
U. 


An' 


Q '70'7 


ioooo 


Q T'QQ 
O. ( OO 


Ar AVlii 


ATA n 

4 ( 4.y 


U.UU ( z 


OOi. 


,i 


i. 


HQ' 

.Uf5, 


o. ( t\) 


1 1 Q'7K 


Q "TQ/I 

o. / o4 


O AVI 


A on c; 

4zy.o 


u.uuiy 


OD. 


c: 
.O 


i. 


1 


o.y4o 


1 Q/l /Ifi 
io44D 


o.y4y 




A Al Q 
44i.o 


U.UUzo 


OO. 


o 
,z 


1 

i. 


Tfi' 

■ ' 'J. 


3.985 


11878 


3.991 


o AVi 


466.8 


U.UU44 


173. 


.4 


1. 


,20. 


A 1 QQ 

4.1cSo 


izUoy 


A 1 QQ 


KjV AVii 


Qon n 

ozy.y 


A AAI O 

U.UUiz 


OA 

yu. 


.o 


A 
U. 


A 1 ■ 
,41 


A '700 


1 1 QTX 


4. ( zy 


D AVi 


A r;'7 A 
40 ( .4 


A ni QA 
U.Uic54 


/I A'7 

4U / . 


.D 


i. 




A "7 a A 
4. / D4 


riQi Q 

yoio 


4. / /U 


Dl AiV 


4Uo.U 


A AAAQ 


1 z, 


o 

.y 


A 
U. 


on" 


4.824 


9810 


4.831 


bl AIV 


441.5 


0.0031 


128. 


,4 


0. 


,48. 


4.941 


9811 


A a AT 

4.947 


Ol AIV 


358.2 


0.0053 


86. 


.8 


1. 


,98. 


o.Uoz 


ii ( o4 


K OQO 

o.uoy 


o AV 


Am o 
4Ui.z 


A nnAA 
U.UU44 


fifi 

DO. 


o 

.z 


o 
z. 


,U/ . 


5.209 


9817 


Ti Ol T 


Dl AIV 


ATO O 

4 ( z.z 


0.0115 


OQT 

zo7. 


.9 


0. 


,92 


5.375 






no ID 




A AAOO 

U.UUzz 


244, 


.9 


0. 


,11 


5.600 


lozooU 


5.606 


iNl AAVi 


332.2 


A AAOQ 

U.UUzo 


62. 


.5 


0. 


,46. 


6.004 






no iiJ 




0.0010 


7o. 


.3 


0. 


1 T 


D.Uio 






no ID 




A AAOK 

U.UUzo 


fi/l 
D4. 


Q 
.O 


A 

u. 


.Ol _ 


D.U4o 


oU4o 


fi oco 


A T YTTT 
Al Aili 


392.0 


A AAfiQ 

U.UUoo 


4i4. 


.3 


0. 


,14 


6.065 






no ID 




0.0014 


55. 


. ( 


0. 


23. 








no ID 




A AAQ 1 

U.UUoi 


c;a 
OU. 


.o 


A 

u. 


,DD. 


D.iUo 


1 QOCOfi 

iozozO 


fi 1 HQ 


iNl AAVi 


OOfi 1 

zzo.i 


U.UUocS 


4U ( . 


1-7 
. ( 


A 

u. 


OT 

,zz 


D.iOO 






no ID 




A AAOA 

U.UUzU 


fio 
Dz. 


y 


A 
U. 


OQ' 
,ZO 


fi 1 K/l 
D.104 


OQI fi 

yoio 


fi 1 QO 

D.icSz 


Ol AiV 


1 QTA T 

io ( 4. r 


A AAQyl 

U.UUo4 


1 fi/1 
iD4. 


.D 


A 

u. 




fi 1 '70 
D.l ( Z 


OQI fi 

yoio 


fi 1 QO 

D.icSz 


Ol AiV 


4yo.o 


A r\A 1 A 

U.U414 


/1 1 n 

4iy . 


,4 


i. 


,oo. 


6.295 






no iD 




0.0075 


332. 


.6 


0. 


,20. 


6.352 


izoibo 


6.360 


re A AiV 


392.0 


0.0138 


425. 


.8 


A 
U. 


,28. 


fi Q'TC 






no ID 




A AAOK 

U.UUzo 


OO. 


c 
.o 


A 
U. 


,O0. 


fi /1 1 Q 

D.4io 






no ID 




A AAOA 

U.UUzU 


OO. 


Q 

.O 


A 

u. 




6.440 


128161 


6.446 


re A AiV 


270.2 


0.0139 


274, 


,1 


0. 


41 


fi AQ(\ 


/ooo 


fi A 

D.4y / 


iVlg Ail 


QQO 7 

ooy. / 


A AA/I Q 
U.UU40 


1 fiO 

ioz, 


7 
. ( 


A 
U. 


1 


6.550 






no ID 




0.0019 


67. 


,4 


0. 


,2o; 


6.569 


128159 


6.575 


Fe XXIV 


260.3 


0.0204 


242. 


.7 


0. 


.74: 


6.590 


7853 


6.580 


Mg XII 


-455.2 


0.0013 


59. 


,5 


0. 


,15: 


6.640 


9691 


6.648 


Si XIII 


361.4 


0.0163 


291. 


,7 


0. 


,42: 


6.729 


7852 


6.738 


Mg XII 


401.2 


0.0100 


54. 


,7 


3. 


.40: 


6.778 


128157 


6.784 


Fe XXIV 


256.7 


0.0269 


289. 


.2 


0. 


,76: 


6.805 


132821 


6.811 


Ni XXVI 


264.5 


0.0198 


416. 


,0 


0. 


34: 


6.842 






no ID 




0.0013 


75. 


,2 


0. 


.ii: 


6.869 






no ID 




0.0090 


271. 


,0 


0. 


,2o: 


6.956 






no ID 




0.0015 


73. 


,4 


0. 


.12: 


7.013 






no ID 




0.0026 


58. 


,3 


0. 


,26: 


7.058 






no ID 




0.0144 


445. 


.5 


0. 


,20: 


7.093 


7855 


7.106 


Mg XII 


557.0 


0.0426 


495. 


,7 


0. 


,6i: 


7.159 


128021 


7.169 


Fe XXIV 


419.1 


0.0556 


392. 


,7 


1. 


17: 


7.223 






no ID 




0.0022 


48. 


,2 


0. 


,25: 


7.242 






no ID 




0.0058 


312. 


.1 


0. 


.11: 


7.265 






no ID 




0.0024 


68. 


,3 


0. 


,18: 





1+0. 


,14 


'-0 


.16 


,+0, 


,07 


'-0 


.06 


,+0, 


,30 


'-0 


.24 


1+0, 


,33 


'-0, 


.28 


.+0. 


,14 


'-0, 


.12 


+0. 


,16 


-0 


.16 


,+0, 


,08 


-0 


.07 


1+0, 


,22 


'-0 


,12 


,+0, 


,14 


'-0, 


.13 


[+0, 


,61 


*-0, 


.37 


,+0, 


,46 


-0 


.37 


1+0, 


,09 


•-0 


.08 


+0, 


,05 


-0 


.05 


;+0, 


,22 


'-0, 


.20 


1+0, 


,09 


'-0, 


.07 


,+0, 


,13 


-0, 


.13 


1+0, 


,02 


^-0, 


.02 


,+0, 


,14 


'-0, 


.13 


;+0, 


,23 


'-0, 


.21 


1+0, 


,02 


■-0, 


.02 


[+0, 


,12 


*-0 


.12 


,+0, 


,04 


-0, 


.04 


;+0, 


,04 


'-0, 


.04 


1+0, 


,03 


'-0, 


.02 


[+0, 


,02 


*-0, 


.02 


;+0, 


,13 


'-0, 


.13 


,+0, 


,13 


-0 


.12 


+0, 


,03 


-0 


,03 


,+0, 


,04 


'-0, 


.04 


1+0, 


,10 


'-0 


.10 


^°o: 


,04 




.04 


■l°o 


,11 




.09 


,+0, 


,03 


■-0 


,03 


•IS: 


,92 




.59 


it°o: 


,03 




.03 


^°o: 


,02 




.02 


+0, 


,07 


-0, 


.07 


't°o: 


,03 




.02 


,+0, 


,08 


■-0, 


.07 


it°o: 


,10 




.10 


't°o: 


,02 




.02 


+0, 


,02 




.02 


,+0: 


,02 


-0, 


.02 


■t°o: 


,12 




.11 


+0, 


,01 


-0, 


.04 


st°o: 


,07 




.06 



lower level upper level ion column^ 



0.014 


ls.2S 


5p.2P° 


8.0111 


X 


1021 


0.412 


1S.2S 


2p.2P° 


i-9tg:i 


X 


1021 


0.029 


ls.2S 


4p.2P° 


9.otl:l 


X 


1021 


0.766 


ls2.1S 


ls.2p.lP* 


'-■°-0.2 


X 


1021 


0.079 


1S.2S 


3p.2P° 


9+0.4 
•^•^-0.3 


X 


1021 


0.416 


1S.2S 


2p.2P'' 


7 4+30 
'•^-2.7 


X 


1021 


0.413 


1S.2S 


2p.2P<' 


r. 4+0-3 


X 


1020 


0.008 


ls.2S 


6p.2P° 


9 9+2.4 
^■^-1.4 


X 


1021 


0.014 


ls.2S 


5p.2P° 


2.9l«:l 


X 


1021 


0.029 


1S.2S 


4p.2P'' 


5.6l}.t 


X 


1021 


0.761 


ls2.1S 


ls.2p.lP* 


4 5+1° 


X 


1020 


0.079 


ls.2S 


3p.2P<' 


Q 1+0-9 
y-J^-0.8 


X 


1020 



0.007 
0.079 

0.024 

0.414 
0.414 

0.002 



0.004 
0.008 



0.079 
0.026 



ls2.2s 
ls.2S 

ls2.2s 

1S.2S 
ls.2S 

ls2.2s 



ls2.2s 
ls.2S 



lS2_7p 



3p.2P° 1. 



,+0.3 
'-0.3 



ls2.5p 7.4: 



+0.8 
0.7 



2p.2P° 
2p.2P° 



2.7 
2.1 



+0.7 
0.6 
+0.1 
0.1 



X 10^ 



X 102 



X lOi 

X 102 



iS2_9p r.stoi X 102 



lS2_8p 
6p.2P° 



6.5 
1.7 



+0.5 
0.5 
+0.4 
0.4 



X 102 
X 102 



0.007 


ls2.2s 


lS2_7p 


6-7lg:! 


X 


1021 


0.014 


1S.2S 


5p.2P'' 


6.9ti1 


X 


1020 


0.748 


ls2.1S 


ls.2p.lP* 


4+0.3 


X 


1019 


0.029 


ls.2S 


4p.2P° 


7 7+2.1 
'•'-1.3 


X 


1021 


0.012 


ls2.2s 


lS2_6p 


Q 7+0.2 
■^•'-0.1 


X 


1021 


0.032 


ls2.2s 


ls2.4p 


7 7+0.5 
'•'-0.4 


X 


1021 



1S.2S 
ls2.2s 



3p.2P'' 
ls2.5p 



4.8: 
2.5 



+0.2 
0.2 
+0.1 
0.1 



1020 



X 10^ 
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Table 2 — Continued 



index 



width>= 



7.283 






no ID 




0.U030 


60.1 


0.27" 


7.373 






no ID 




0.0051 


47.7 


0.64" 


7.387 






no ID 




0.0029 


65.1 


0.22" 


7.400 






no ID 




0.002U 


61.6 


0.15" 


7.415 






no ID 




0.0029 


61.9 


0.23" 


7.432 






no ID 




0.0026 


54.2 


0.24" 


7.464 


126341 


7.472 


Fe XXIII 


321.5 


0.0272 


348.1 


0.48" 


7.490 






no ID 




0.0036 


58.0 


0.31" 


7.513 






no ID 




0.0023 


52.0 


0.21" 


7.533 






no ID 




0.0034 


59.9 


0.27" 


7.555 






no ID 




0.0020 


59.7 


0.15" 


7.575 






no ID 




0.0042 


53.8 


0.38" 


7.617 






no ID 




0.0022 


43.5 


0.23" 


7.664 






no ID 




0.0067 


278.2 


0.10^ 


7.700 


126899 


7.733 


Fe XXIII 


1285.7 


0.0027 


48.3 


0.25" 


7.722 


120899 


7.733 


Fe XXIII 


427.3 


0.0018 


46.4 


0.17" 


7.743 


126899 


7.733 


Fe XXIII 


-387.5 


0.0030 


72.0 


0.18" 


7.850 






no ID 




0.0056 


48.5 


0.55" 


7.863 






no ID 




0.0037 


53.6 


0.29' 


7.875 






no ID 




0.0037 


57.2 


0.27' 


7.890 






no ID 




0.0028 


51.8 


0.22 


7.904 






no ID 




0.0035 


59.5 


0.24 


7.946 


127946 


7.983 


le XAlV 


1396.9 


0.0019 


51.1 


0.14' 


7.980 


127946 


7.983 


Fe XXIV 


112.8 


f\ noon 

0.0829 


409.0 


1.24 


8.081 


120327 


8.090 


Fe XXII 


348.9 


0.0043 


52.5 


0.33" 


8.296 


125709 


8.303 


Fe XXIII 


249.5 


0.0431 


311.2 


0.67" 


8.393 


7883 


8.421 


Mg XII 


1000.8 


0.0042 


50.3 


0.29" 


8.409 


7883 


8.421 


Mg XII 


428.1 


0.0704 


285.7 


1.28' 


8.705 


119538 


8.715 


Fe XXII 


344.6 


0.0061 


54.4 


0.36^ 


8.721 


119538 


8.715 


Fe XXII 


-206.4 


0.0057 


57.7 


o.3i; 


8.963 


119527 


8.977 


Fe XXII 


478.6 


0.0077 


41.5 


o.6i: 


9.048 


132824 


9.061 


Ni XXVI 


431.0 


0.0602 


244.6 


0.88: 


9.092 


132807 


9.105 


Ni XXVI 


428.9 


0.0409 


180.2 


0.76: 


9.158 


7744 


9.169 


Mg XI 


353.8 


0.0038 


36.1 


0.28: 


9.336 


132741 


9.340 


Ni XXV 


118.9 


0.0414 


304.6 


0.42: 


9.373 


132694 


9.390 


Ni XXV 


534.5 


0.0289 


164.9 


0.48: 


9.469 


6241 
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0.2 



X 10^ 
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X 
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2 5+1-8 


X 
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X 


1020 
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X 
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X 
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X 
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Table 2 — Continued 



index 



width'^ 



10.490 


124935 


10.505 


Fe XXIII 


443.3 


0.0233 


66.5 


0.77' 


10.570 


127980 


10.619 


Fe XXIV 


1390.7 


0.0309 


479.7 


0.12' 


10.610 


127980 


10.619 


Ic aXIV 


254.5 


f\ o o o o 

0.2323 


490.0 


1.30' 


10.650 


127947 


10.663 


be AAlV 


366.2 


0.3370 


716.3 


1.13' 


10.780 


124988 


10.785 


Fe XXIII 
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0.0346 


478.1 


0.12" 


10.890 
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10.903 
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352.6 


0.04U2 
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Fe XXIII 


1372.3 
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41.0 


0.31' 
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0.74" 
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11.018 


l*e XXI 11 
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0.1051 


221.2 


0.99' 


11.150 






no ID 




0.0350 


382.1 


0.15" 


11.280 


125155 


11.299 


Fe XXIII 


508.0 


0.0289 


253.5 


0.15" 


11.310 
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11.299 


Fe XXIII 


-289.1 


0.0224 


170.2 


0.18' 


11.350 


124975 


11.337 


Ic AAlii 


-348.9 


0.0504 
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0.16' 


11.420 
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11.442 
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585.8 


0.1575 
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0.52' 
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11.457 


Fe XXIII 


-337.4 


0.0621 


253.4 


0.34' 


11.540 






no ID 
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437.0 


0.11' 


11.760 


119502 


11.769 


Fe XXII 
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0.0710 


192.6 


0.51' 


11.910 


120311 
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Fe XXII 
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0.0546 
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0.50' 
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Table 2 — Continued 





index 




ion 


vo° 


ew'' 


width-^ 




ff 


lower level 


upper level 


ion column^ 


14.550 
14.610 
14.770 
14.820 


4675 
4634 
4634 


14.645 
14.832 
14.832 


no ID 
O VIII 
O VIII 
O VIII 


724.9 
1255.2 
238.9 


0.0577 
0.1763 
0.0501 
0.1715 


26.3 
189.3 

26.7 
454.9 


7 7-1+9.99 
'■'-'--4.86 
75+0.30 

3 n2+«-"« 
•J-U^_1.90 


0.008 
0.014 
0.014 


ls.2S 
ls.2S 
ls.2S 


6p.2P° 
5p.2P° 
5p.2P° 


1.8+°-^ X 1020 

4.1+9-9 ^ ;^q20 

3.5+i° X lOis 



^Wavelength in A 

''Wavelength taJcen from the xstar database. 

'^Velocity is in units of km see" ^ 

''Equivalent width in units of eV. 

''ine center optical depth. 

^Oscillator strength 

^Column density is in units of cm"^ 
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The range between 1-2.5 A is dominated by the H-like Lq and He-like allowed n =1 - 2 lines 
from the elements Fe (1.77, 1.85 A ), Mn (1.92, 2.01 A ) and Cr (2.08, 2.18 A ). Wavelengths 
shortward of the Fe XXVI La line at 1.77 A contain the K lines of elements heavier than Fe, in 
addition to the higher series lines of Fe. Although there are marginal detections of some of these 
lines, there are too few counts in this region to strongly constrain the abundances of the heavier 
elements Co, Cu, Zn. Ni is an exception to this, since it also has lines from Li-like Ni XXVI at 
longer wavelength. It is important to emphasize that, although Lya lines are not clearly detected 
for Co, Cu and Zn, we cannot set meaningful upper limits to their strength because our model 
atoms for these elements do not include n =2 - 3 lines from the Li-like and lower ionization stages. 

The 2 - 3 A region includes the He-like resonance line of Mn XXIV, the K lines from H- and 
He-like Cr (2.088 and 2.179 A ), He-hke Sc (2.877 A ), H-hke Ti (2.491 A ). The apparent absence 
of lines from V, H-like Sc and He-like Ti provides relatively secure upper limits on the abundances 
of these elements. We point out that the existence of evaluated wavelengths for these lines from 
NIST provides added validity to this conclusion. In this range are also n =1 - 3 lines from Ca, 
which are stronger than the n =1 - 2 lines from Sc, Ti, Cr and Mn. In this range are also n =1 - 
3 lines from Ca, Ar, S, Si. For Si XIV Rydberg series lines are detected up to n = 6. 

At 4.18 A we marginally detect the La line from H-like CI XVII. The He- like line from this 
element is not detected, but we note the absence of evaluated wavelength for this line. Similar 
comments apply to the H-like line from P, near 5.38 A. The ratio of the He and H-like lines from S 
near 4.72 and 5.04 A follows the same behavior as for Ca and Ar. The 6 -7A region is dominated 
by lines from Si XIII and XIV and Li-like Fe and Ni; for Fe these are detected up to ls^2s — Is^lOp. 
At 10.67 and 10.88A are the lines from Fe XXII which provide the density diagnostics. These lines 
together extend the range of ionization of observed species to include ions which are indicative of 
lower ionization parameter gas than the hydrogenic and helium-like ions which predominate. 



Table |2] includes 175 hues of which 100 have IDs. The table of [Miller et"aL] pOOSl ) includes 



102 lines, of which 15 have no IDs. Our table has 44 lines which are not in the Miller et al. (20081 
table, although we note that they used a 5a criterion for line detection, while ours is 3a. Of the 
unidentified lines in 



Miller et al. 



(20081 we propose IDs for 5.6 A Ni XXVI 2s - 7p and 9.372 A 



Ni XXV 2s^ - 2s3p. We have no IDs for the unidentified hues from [Miller et al.| ([2008[) at AA 
6.25, 7.0555, 7.0851, 9.9509 A. These latter 3 are crudely consistent with lines arising from the 2p 
excited levels of Ni XXVI (rest wavelengths 7.048, 7.0950 and 9.6383). But the results of warmabs 
suggest optical depths for these lines which are smaller than the resonance lines by factors > 10^ 
assuming a density of 10^^ cm~^ and no trapping of line radiation. All the 75 lines from our notch 
model for which we have no identifications, have small equivalent widths, within a factor of 2 of 
our cutoff of 3.5 xlO"*^. We speculate that some of them could be due to the notch algorithm 
attempting to correct for discrepancies between the model continuum and the data, perhaps due to 
bound-free continuum absorption. In addition, visual inspection of the spectrum shows lines which 
do not appear in our table. Many of these have possible IDs as lines from excited levels, but none 
appears in the warmabs models at sufficient strength to be included in our table. Other lines in 
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this category include: 2.60 A possibly Ti XXI He-like 1-2, 2.15 A possibly Sc XX He-like 1-3, 
features at 3.61, 5.375, 5.62, 5.65 A with no obvious ID, and 5.72 A possibly Al XIII 1-4. Features 
at 6.22 and 6.24 A could be due to Si XIII ls2s - 2s2p. Discernible features at AA 6.7, 6.87, 7.06, 
9.24, 9.78, 9.96, 10.1, 10.15 and 11.15 A have no obvious IDs. Many of these correspond to lines 
in our list which are from excited levels, or from ion stages which are too low to coexist with the 
dominant ions in the WARMABS models. We also note an emission line at 13.38 A which coincides 



with the A 13.387 A 2s - 3p transition in Ti XXII listed in Shirai et al. (2000). If this is correct, 
it is hard understand why this line should appear in emission when other analogous lines from Fe 
XXIV and Ni XXVI, for instance, appear in absorption. Other interesting lines not identified by 
Miller et al.l (|2008|) are the Fe XXIII hues between 11.28 A and 11.47 A, since they arise from the 



metastable 2s2p P level. Their presence corroborates the density estimate from the Fe XXII lines. 



The Doppler shift of the lines relative to the lab wavelength are shown graphically in figure 15 



In this figure, the color corresponds to the element, and the size of the dot corresponds to the line 
optical depth. This shows that the lines cluster in a range around 400 km s^^ from zero offset (here 
and in what follows we quote blueshifted velocities as positive, and conversely). This corresponds 
to ~2 - 10 mA for most lines. Many of the laboratory wavelengths are uncertain at this level. 
Given this, the most notable aspect of the line shifts is the large velocity offset of the Fe XXVI La 
and the Ni XXVII lines. These are among the strongest lines in the spectrum, and are blueshifted 
by ~1300 km/s, which is significantly greater than any other lines in the spectrum. This, together 
with the fact that these are the two highest-ionization lines in the spectrum, suggests that these 
lines are partly formed in a separate component of the flow. If so, this component would have 
higher ionization, and higher velocity, than the component responsible for the rest of the lines. On 
the other hand, we consider the possibility that this is an artifact of shortcomings in the HETG 
calibration when applied to fitting absorption of the steeply sloping continuum under these lines. 
In what follows we will examine alternative explanations for this in our fits to the spectrum. 

The line widths in Table [2] range from 1-10 mA. We have searched and found no significant 
correlation between this quantity and simple quantities characterizing the parent ion, such as the 
ionization potential or isoelectronic sequence. There is a weak tendency for the largest widths 
to be associated with lower ionization potential ions. An example is the Lq line from Si XIV. 
However, there are also narrower lines from ions with comparable ionization potential, so this does 
not constitute a statistically significant correlation. 



2.2. Ion Column Densities 



Ion column densities can be derived from the line identifications in Table [2] assuming that the 
equivalent widths lie on the linear part of the curve of growth. That is, we calculate 
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Line Velocity Shifts 
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Fig. 15. — Doppler velocities obtained by comparing measured notch wavelengths with lab wave- 
lengths from XSTAR database. Element color coding is given in the legend. The dot size is propor- 
tional to the line optical depth. 



N- - '^''"'^ (I) 

where Avo is the total Doppler width in frequency units, including both thermal and turbulent 
contributions. is the line depth measured from the notch fit to the spectrum, and since it is 
based on Gaussian line fits its accuracy diminishes when values become large. In Table |2] we list the 
value of the ratio iVjon/^eiement) where VeZement is the elemental solar ( Grevesse et al.|[l996 Allen 



1973) abundance relative to hydrogen. Nion/yeiement, is the equivalent hydrogen column density 
implied by a given line if its ion fraction were unity, if the elemental abundance were solar, and if 
it lies on the linear part of the curve of growth. In figure [16] we plot this quantity. The horizontal 
axis is related to the atomic number of the parent element by Zion = Z^iement -|- 1 — 0.1 * (ionstage), 
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where {ionstage)=l for hydrogenic, 2 for He-like, etc. Error bars are plotted and, in most cases, 
are small compared with the the interval between points. If the assumptions of unit ionization 
fraction, cosmic abundances, and linear curve of growth were correct, then all the points in this 
diagram would lie along one horizontal line. Since most elements have lines from more than one 
ion, the assumption of unit ionization fraction cannot be correct; for these the true column should 
be the sum of the columns for various ions if the other assumptions were correct. However, this 
neglects the possible contributions from ions which do not produce observed lines, such as fully 
stripped species, which are likely to be important for the lower-Z elements. 



It is clear from figure 16 that there is a greater dispersion of column densities within elements 
than can be accounted for by ionization effects. This will be discussed in greater detail in the 
next section, and suggests the limitation of our second assumption, that of the linear curve of 



growth. Figure 16 also shows departures from solar abundance ratios, in the sense of apparently 
enhanced abundances for most elements between Sc and Mn, relative to elements with Z < 16. This 
conclusion is dependent quantitatively on the excitation and ionization conditions in the absorber, 
and in the following section we will attempt to quantify the abundances for a variety of assumptions 
about the state of the gas. 

We expect qualitatively that the lower-Z elements will be more highly ionized than the higher-Z 
elements, for most plausible ionization mechanisms. This would predict that the apparent elemental 
abundances would be systematically greater for the low-Z elements than for the high-Z elements in 
figure [T6j since the low-Z elements would have a greater fraction of their ions in the unobservable 
fully stripped stage. This is the opposite behavior to what we observe, and so reinforces the 
conclusion that the elements with Z > 16 have enhanced abundances relative to those with lower 
Z. 



2.3. Curve of Growth 

The validity of the procedure used to derive the columns in figure [T6| depends on the assumption 
that the lines lie on the linear part of the curve of growth. If this is correct, then the line equivalent 
widths should be proportional to the transition oscillator strengths, and comparison of various lines 



from the same ion should show this proportionality. In figure 17 we test this procedure using the 
Lyman series lines from H-like ions of Ne, Mg, Al, Si, S, Ar, Ca, Cr, and Mn. We also include the 
2s — np lines from Li-like Fe XXIV and Ni XXVI. Solid lines are linear regression fits to data with 
errors less than 0.1. Also shown on this figure are diagonal lines (dashed) corresponding to the 
proportionality expected for linear curve of growth. This shows that, although some ions appear 
to follow the linear trend, the strongest lines in particular grow more slowly than linearly. This is 
particularly apparent for the lines of Ne X, Fe XXIV and Ni XXVI. Each ion has at least 5 lines, 
and the trend is apparent across the line strengths. This shows that the simple analysis provided 
in figure 16 and Table [2] are likely not adequate for the purposes of inferring abundances and the 



lines may be saturated. On the other hand, weaker lines such as those of Si XIV do apparently 
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follow the linear trend. 



Another illustration of the effects of saturation is shown in figure 18 This shows the ratios 
of He-like to H-like ion abundances inferred from the line equivalent widths and identifications in 
Table [2j These ratios are independent of element abundance, owing to the fact that each ratio 
is taken between ions from the same element. This shows that the ratios are all in the range 
between 10"^ and 10+^ for 12 < Z < 26, but that there is no systematic trend with Z. We might 
expect the ratio to increase with Z, as higher Z elements would be generally less highly ionized. 



for plausible ionization mechanisms. Also shown in figure 18 are the contours traced by an xstar 
photoionization model described in more detail in the Appendix. These are labeled by the value of 
log(^), where ^ = L/{nE?) is the ionization parameter as defined by Tarter et al. (1969); L is the 



source energy luminosity integrated from 1 - 1000 Ry, n is the gas number density, and R is the 
distance from the continuum source to the absorber. These show that a given value of ionization 
parameter predicts that the He/H abundance ratio should increase between adjacent elements by 
a factor ~5. The figure clearly shows that a single ionization parameter cannot account for the 
ratios displayed by all the elements. This could indicate the existence of a broad range of ionization 
parameters in the source, spanning values indicated by this figure. On the other hand, it could 
be associated with radiative transfer effects, such as saturation, which make ion fractions inferred 
from a linear curve of growth unreliable. 

Possible explanations for the departures from the linear curve of growth include the influence 
of saturation which causes the curve to flatten when the lines become optically thick in the Doppler 
core. Other possibilities include filling-in of the lines by an additional continuum emission com- 
ponent which is not seen in transmission through the warm absorber, and also radiation transfer 
effects in the absorber itself. The latter includes forward scattering, which may also depend on the 
relative size of the continuum source and the absorber. Also, thermal emission can fill in the lines. 
In the following subsection we discuss these in turn. 



2.4. Radiation Transfer and Curve of Grovi^th 

The standard curve of growth for a resonance line can be written: 



EW = de(l- e-^'^^A (2) 



r{e) = ^f^^^Y^N-^<p{-^) (3) 

rrieC Vturb ^^turb 

where e is the photon energy, / is the oscillator strength, Xij is the ion fraction, Yj is the element 
abundance, N is the total column density, A is the line wavelength, VTurb is the turbulent velocity 
(including the thermal ion speed), and (j) is the profile function, which includes both a Doppler core 
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and damping wings. This is shown in figure 19 for various choices of the velocity characterizing 



the Doppler broadening vturb-, and for a natural width corresponding to that for the Si XIV La 
line. This shows that, for a given equivalent width, the effects of saturation are more likely to be 
important when the Doppler broadening is smallest. That is, the equivalent width where the curve 
flattens is approximately proportional to the Doppler width. As we will show, the widths of many 
lines from GRO J1655-40 are not constrained from below by the HETG spectrum. Thus, a possible 
explanation for the slower than linear curve of growth is that the line Doppler widths are small 



enough such that the strongest lines in figure 17 are affected by saturation. 



Filling in by continuum from a separate source which is not absorbed by the warm absorber 
would preferentially affect lines with larger optical depths. This could qualitatively explain the 
flattening of the curve of growth which is observed. However, the filling in would also have the 
wavelength dependence of the additional continuum component. Although this is unknown, the 
simplest assumption would be that it is the same as the primary continuum. We have tested this 
possibility quantitatively using direct fitting, and will discuss this further in the following section. 

The standard curve of growth in equation ^ assumes the simplest possible geometry and 
microphysics affecting the line: (i)The continuum source has negligible size; (ii) The absorber 
exists only in a narrow region along the line of sight to the continuum source; (iii) The excitation 
and deexcitation of the transition responsible for the line is affected only by the radiative excitation 
and spontaneous decay connecting the two atomic levels. In addition, the gas is assumed to have a 
velocity field which can be characterized by a Maxwellian distribution with a well-defined thermal 
or turbulent velocity. 

If the continuum source has finite size, or if the warm absorber exists outside the line of sight 
to the source, then the line profiles will be affected by photons which scatter into our line of sight. 
Although there is no simple general expression for the intensity observed from an extended scatter- 
ing atmosphere in this case, it is straightforward to show that in the limit of small optical depth, 
the scattered emission contribution scales proportional to [D/R)'^, where D is the characteristic 
size of the continuum source and R the characteristic distance from the continuum source to the 
absorber. This contribution is independent of optical depth in this limit, and so will not affect the 
shape of the curve of growth at small r. 

Departures from assumption (iii), concerning the population kinetics, can be expected if the 
radiative decay of the upper level is suppressed or if the upper level can be populated by another 
process such as recombination or collisional excitation. Some of these possibilities were discussed 



by Masai Sz Ishida (2004 1 . If the upper level can decay via alternate channels, such as branching 
to other levels radiatively, collisional deexcitation, or autoionization, the basic absorption profile 
properties will be unchanged (although these processes can affect the damping parameter). Sup- 
pression of the upper level decay by collisional deexcitation could occur if the product of line optical 
depth and electron density, ng, is large, but this requires optical depth ~ lO^^/rie for iron. This 
would affect our limits on possible filling-in of the line by collisional or recombination radiation 
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which we present in the next section. Populating the upper level by collisions or recombination 
requires suitable temperature and ionization conditions. Efficient population by coUisional exci- 
tation requires temperatures greater than can be accounted for by photoionization heating and 
radiative cooling, and this in turn will affect the ionization balance. It also implies the existence 
of an additional heating mechanism. We have tested these possibilities quantitatively using direct 
fitting, and discuss the results below. 



3. Direct Fitting 

Many of the issues discussed in the previous section can be tested using direct fitting. These 
include: detector resolution, counting statistics, and scattered and thermal (either collisional or 
recombination) emission. In order to do so, we use the WARMABS analytic model which interfaces 
with the XSPEC spectral fitting package. WARMABS makes use of a stored, precomputed table of 
level populations for a family of photoionization models. It uses these to calculate a synthetic 
spectrum 'on the fly' within xspec. The advantage of this is that it calculates the synthetic spec- 
trum automatically for the energy grid of the data such that it matches the energy resolution of 
the instrument. Interpolation can introduce significant numerical errors when applied to absorp- 
tion spectra, since the absorption coefficient can change rapidly over a narrow range in energy. 
Warmabs calculates all absorption lines using a Voigt function including damping due to both 
radiative and Auger decays where applicable. Line profiles are calculated using sub-gridding, on an 
energy scale which is a fraction of a Doppler width, and the opacity and transmittivity then mapped 
onto the detector grid. Bound-free absorption is also included. WARMABS uses the full database 



and computational routines from the xstar package (Kallman and Bautista 2001 Bautista and 
Kallmaiil|2001 1 , and differs from xstar in that the level populations are pre-calculated rather than 
calculated simultaneously with the spectrum. The level populations are calculated using a full 
'collisional-radiative' calculation, albeit with relatively simple model ions in most cases. Thus they 
include the effects of upper level depletion due to thermalization and photoionization automatically. 
We have extended the xstar database to include all the trace elements seen in the GRO J1655-40 
spectrum. A description is provided in the Appendix. 

Warmabs does not calculate the flux transmitted by a fully self-consistent slab of gas, as xstar 
does. Rather, it calculates the opacity from the illuminated face of such a slab, and then calculates 
optical depth and transmitted flux analytically by assuming the opacity is uniform throughout 
the slab. A real slab will shield its interior, which will then have lower ionization and greater 
opacity than a slab whose opacity is assumed to be uniform. Thus, the results of warmabs fltting 
will be absorbing columns which may be greater than would be produced by a real slab. The 
importance of this effect depends sensitively on the assumed slab structure, i.e. its geometrical 
thickness, and so its importance cannot be estimated in general. Also, this effect is negligible at 
high ionization parameters. As we will show, our flts require large ionization parameters. In this 
case, the advantages of the energy resolution provided by warmabs outweigh its limitations. 
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3.1. Model 4: Narrow Lines 

Our basic fit assumes that the absorption is provided by a single component of photoionized 



gas. We adopt the ionizing continuum used by Miller et al. (20081 which consists of a power law 



plus disk blackbody, as described in the previous section. The power law photon index is 3.54 
and the disk inner temperature is kT=1.35 keV. In our fit we account for the curve of growth 
by adopting a small turbulent velocity, 50 km s~^, so that the ratios of these lines are on the 
saturated part of the curve of growth. We then search for the single ionization parameter which 
most nearly accounts for all the lines in the spectrum. We do this using warmabs and the xspec 
package and stepping through ionization parameter in intervals Alog(^)=0.1. Warmabs accept the 
turbulent velocity as an input, but the actual line width is calculated including both this turbulent 
velocity and the thermal ion velocity corresponding to the equilibrium temperature. As shown in 



the Appendix figures 23 - 26 in order to simultaneously produce the ions of iron ranging from 
B-like (Fe XXII) to H-like, an ionization parameter in the range 3<log(^)<4 is needed. We find the 
best fit at log(^)=4. 0^0.1, log(N)=23.8 and a blueshift for the absorber of 375 km s~^. This can be 
compared with the radial velocity of the system of 141^1 km and the orbital semi-amplitude 



of the secondary star of 215.5 ^2.4 km s ^ (Shabhaz et al. 2002 1 . This suggests that the outflow 



is not moving fast when compared with the maximum velocities characterizing the orbital motion 
in the system. We have also performed equivalent fits for this set of assumptions using the ISIS 



fitting package (Houck & Denicola 20001, and have verified that the results are independent of 



which fitting package is used. 

We refer to this single component warmabs with Vt^rb = 50 km s""*^ as model 4. We assume 
a gas density of 10^^ cm~^ in calculating the level populations used by warmabs. We have not 
extensively tested our results at lower densities, and we point out that it is likely that densities 
as low as 10^^'^ cm~'^ can produce the Fe XXII lines. We employ the same ionizing continuum as 
derived in the model 1 fit when calculating the populations and gas temperature using xstar. In 
our fitting for model 4 and subsequent models we allow the normalizations of the two continuum 
components to vary. The best fit flux is 1.8 xlO~® erg cm~^ s~^ in the 2-10 keV band. In Figures 

[HI 

we show the observed count rate spectrum (black) together with all the models discussed in 



this section as red curves, labeled according to the model number. Strong lines are marked in blue 
along with parent ion. 



As discussed above, and shown in figure 15 the majority of lines have centroid wavelengths 



which are consistent with a single Doppler blueshift, approximately 400 km with respect to 
their laboratory values. The notable exception is the Fe XXVI La line near 1.77 A which has a 
centroid wavelength corresponding to a Doppler blueshift of ~1300 km s~^. In model 4, we adopt 
the hypothesis that this line arises in a separate velocity component of the flow, which must have 
a much greater ionization parameter so that it does not show up in the other lines. Since the 
XSTAR grid of model level populations does not extend to values where Fe XXVI is the only ion 
of significant abundance, we do this by adding a single Gaussian component for the high velocity 
part of this line. A consequence of this is that the 375 km s^^ component of the flow, which we 
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model with warmabs and an ionization parameter log(^)=4, accurately fits to the red part of the 
Fe XXVI La profile, while the Gaussian accounts for the rest. This is shown in figure 20 which 
shows the region surrounding the line. The data is shown as the black bars, and model 4 is shown 
as the dark blue points. The green points show the continuum + Fe XXVI La Gaussian model 
alone. The equivalent width of the Fe XXVI Gaussian is 29.6 eV. The for this model is 34211 
for 8187 degrees of freedom. These parameters are summarized in Table [T| 

An alternative possibility is that the Doppler blueshift of the Fe XXVI line is affected by 
instrumental effects and shortcomings in the available HETG response function. This could be due 
to the fact that the continuum count spectrum is steeply decreasing in the line region, as shown 
in figure 20 Internal scattering could cause photons from the continuum adjacent to the longer 
wavelength wing of the line to scatter into the line core, and this effect would be stronger on the 
long wavelength side than on the short wavelength side. The width of the line spread function 
is comparable to the line blueshift. In addition, the wavelength calibration may be affected by 
the use of continuous clocking mode. On the other hand, other lines in the 1.5 - 2 A region, 
although at slightly longer wavelength, do not show this effect, as exemplified by the Fe XXV 
Is^ — ls2p line. Furthermore, we have confirmed that the response matrix we use reproduces the 
line spread function obtained from the Chandra calibration database So, although there is no 
obvious shortcoming in the response matrix which would explain the discrepancy between the Fe 
XXVI La wavelength shift and those of other lines, we will examine both hypotheses: that the 
high velocity part of the Fe XXVI line is associated with a separate kinematic component (as in 



model 4), and that it is not. The latter case also corresponds to the assumptions of Miller et al. 



(2008 1 . The primary conclusions of this paper, with respect to dynamics, abundances, turbulence, 
and geometry, are not dependent on this assumption. 

The model 4 fit accurately reproduces the strength and shape of the Fe XXVI and Fe XXV 
lines, although we predict a slightly stronger red wing on the Fe XXV than is observed. This is the 
contribution of Fe Ka from lower ion stages of iron, Fe XXIV and XXIII. The strengths of these 
lines provides a lower bound on the ionization of iron. We also slightly overestimate the strength 
of the Mn XXV La line, relative to the Mn He-like Ka line. This may be due to an error in the 
ionization balance for this element; as shown in figure 18 there are apparent departures from a 
simple monotonic trend in the H/He ratio with nuclear charge. The region longward of 3A includes 
the lines from Ca (3.02, 3.18 A) and Ar (3.73, 3.94 A). Our single component ionization balance is 
too low for both these elements, adequately accounting for the He- like lines but under-predicting 
the H-like lines. We also point out that our assumed line width of 50 km s^^ adequately fits the 
profiles of essentially all the lines shortward of 6 A. 

Model 4 is based on a physical model for the absorber in GRO J1655-40, and it represents 
the best fit that we can obtain to the spectrum using a single warmabs component. As shown in 
figure T]|T4 this accounts for the depths and positions of essentially all the lines in the spectrum. 



^http:/ /cxc. harvard.edu/caldb/calibration/gratings. html 
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Discrepancies between the model and the data fah in three categories: errors in the continuum, 
errors in hne widths, and errors in hne strengths (including lines which are missing). In what 
follows we will explore these in turn. As we will show, the best fits we obtain to a physical 
model have jv ~ 2.9. This can be compared with the phenomenological notch fit, model 3, 
which has reduced jv ~ 2.2. The difference likely represents the ability of the phenomenological 
notch fit to account for features which may not be physical lines but rather part of the continuum 
or consequences of calibration uncertainties. This, we think, represents the ultimate limit in our 
ability to fit the spectrum. In what follows, we proceed and attempt to interpret the values from 
physical models when compared with each other, but without relying on standard interpretations of 
these values and their relation to probability of random occurrence, etc. That is, we acknowledge 
that our fits are not acceptable based on these standard arguments about x^, but nevertheless 
interpret confidence intervals on the fitting parameters by evaluating Ax^ as if they were. 

Model 4 is designed to fit to the curves of growth of lines from Fe XXIV and Ni XXVI by 
having a relatively small turbulent velocity, so that the lines from these ions are at least partially 
saturated. This works for these ions, as is apparent from figure |6] and it also is consistent with 
the observed widths of the lines which are not resolved. However, it does not account for the 
observed widths of the broader lines, such as Si XIII and Si XIV, and also for the 2s - 3p lines of 



Fe XXIV near 10.6 A. [Miller et al^ ( |2008[ ) adopt a hne width of 300 km s~^ which is a better fit to 
the observed widths for many lines. We conclude that saturation with a single profile component, 
although it accounts for the curves of growth of the lines of Fe XXIV, is not consistent with the 
curves of growth of lines from Si XIV, nor is it consistent with observed widths of some of the strong 
lines. A possible explanation is that the lines consist of multiple narrow components closely spaced 
in velocity so they appear blended together in the HETG, and where the components furthest from 
line center are unsaturated and so appear only in the low members of the Li-like 2s - np series. UV 



warm absorber lines in AGN appear to follow this behavior (Gabel et al. 2005). 



3.2. Model 5: Two components plus broad lines 

Another way to fit the Fe XXVI La line is to add a higher ionization parameter warmabs 
component to the fit. We do this in model 5, which has a warmabs component with the same 
ionization parameter as model 4, but with the addition of a second component at an ionization 
parameter which can make the Fe XXVI line sufficiently strong without over-producing the lower 
ionization lines. This model also has a turbulent velocity of 200 km s~^. We vary the column 
densities and element abundances of both components in order to find the best fit. In doing so, 
we force both components to have the same bulk outflow velocity, the same turbulent velocity, 
and the same elemental abundances. Thus, the two components are allowed to differ only in their 
ionization parameters and column densities. This produces a fit with 'x^ jv = 36281/8189. This fit 
is also shown in figures [T]- [14] In figure [20[ we show a blowup of the iron K region, comparing this 
model with model 4. Model 4 is decomposed into the part accounted for by the Gaussian (green) 
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and the total (blue). Model 5 is shown as the red points. This shows that in spite of the nearly 
equivalent x^, model 5 does less well than model 4 in fitting the iron lines; it generally over- predicts 
their strengths. This is particularly true in the case of the Fe XXVI La line. Since this model has 
the same 375 km s~^ outflow velocity as model 4, the best fit accounts for the blue edge of the 
Fe XXVI La by over-predicting the line near line center and the red edge. Thus, the conclusions 
from this model differ from model 4 (and subsequent models) in that they do not depend on the 
existence of two kinematic components. 

In other ways, the comparison between models 5 and 4 reflects the fact that model 5 has on 
average a higher ionization parameter. As a result, model 5 tends to over-predict the ratio of H-like 
to He-like lines, compared with both model 4 and with the observation. Also, model 5 has a larger 
turbulent velocity, and therefore predicts a steeper curve of growth for most lines. This fact is 
apparent from the Fe XXIV lines in the 6 - 7 A region. On the other hand, model 5 fits better than 
model 4 to lines which appear to follow the linear curve of growth, such as the L/3/La lines of Si 
XXIV, and also to lines which are detectably broadened, such as the 10.6 A doublet of Fe XXIV. 

3.3. Model 6: Partial Covering 

An additional possible reason for the apparent saturation of the curves of growth could be due 
to systematic effects or calibration uncertainties associated with the Chandra telescope, grating 
and detector. This could lead to apparent scattering of photons in the continuum into the cores of 
absorption lines which is not accounted for by the detector response matrix, thus preventing the 
residual flux in the lines from going to zero. This effect cannot be evaluated accurately without 
calibration data which includes narrow absorption line features, but we can get an indication of its 
plausibility by fitting to a model which includes some 'leakage' of photons into lines from adjacent 
continuum regions. We do this by fitting to a partial covering model in which the WARMABS 
component is partially diluted, i.e. model = {{1 — C) + C x warmabs) x continuum where 1 — C 
is the fraction of scattered continuum at each energy. We also adopt a turbulent velocity of 200 
km s~^ for this model, since we do not need to account for saturated curves of growth; the dilution 
of the WARMABS lines tends to flatten the curves of growth. We show this as model 6 in Table [T] 
and in figures |l]-[l4| This model has xV'^ = 24561/8184, and so is the best of the fits using the 
WARMABS model. It also employs a turbulent velocity which is sufficient to account for the widths 
of resolved lines while at the same time fitting to the curves of growth for Fe XXIV and Ni XXVI. 

The best-fit value of C is 0.37. We can interpret this as being due to true partial covering in 
the object, i.e. if the continuum source is more extended than the absorber, or as being due to 
instrumental scattering which is not accounted for by the response matrix we used. The latter can 
be further subdivided into shortcomings in the response in accounting for the line spread function 
(LSF) for the HETG which occur in the core region of the LSF, and shortcomings which occur in 
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the wings. The core region of the LSF [5 can be approxmiately represented by a Gaussian with 
full-width-half-maximum (FWHM) of ~1300 km s"^ at the wavelength of the iron line, 1.77 A. So 
two intrinsically narrow lines of equal intensity will fill the region between them to 30% of their 
peak intensities if they are separated by ~ 900 km s~^. This is comparable to the typical line widths 
we find in our fits. However, as discussed previously, the response matrix we use does accurately 
reproduce the LSF in the core region, so it would require a large error in the calibration files, the 
LSF and the resulting response matrix, in order to account for the filling in of the lines in this way. 

Another possibility is that there exists significant scattering in the line wings which is not 
accounted for by the calibration. This is difficult to evaluate quantitatively, except to point out 
that the uncertainty in the LSF at the extremes of the wings appears to be at most ~2 - 4 % of the 
maximum value. In order to make up the covering fractions we require, the integrated area under 
the wings would have to be ~30 - 40 % of the area in the core of the line. The wings would have 
to extend to ~10 times the a of the core, or ~7600 km s~^. We cannot evaluate the likelihood of 
this possibility reliably, and cannot conclusively rule it out, although such a large departure from 
the calibration would be surprising. This observation differs from typical observations in its use of 
continuous clocking mode. This prevents the use of detector regions adjacent to the readout strip 
for background subtraction, but the counting rate in this case is high so that background should 
not be important at the levels considered here. The fact that the response matrix agrees with the 
LSF in the core of the line is further evidence that the use of continuous clocking is not the source 
of the apparent partial covering. 

We conclude that true partial covering in the source is more likely than instrumental effects, 
and thus plays a role in determining the line curves of growth and the overall quality of the fit. 
Owing to its superior we adopt this model as the one which most nearly accounts for the observed 
spectrum and consider the physical assumptions to be most nearly correct for the GRO J1655-40 
outflow. We emphasize that the partial covering does not account for the high velocity component 
in Fe XXVI La, and model 6 includes the same high velocity Gaussian contribution to this line as 
model 4. 

3.4. Other Models 

We also have examined the possible influence of thermal emission filling in lines on the one- 
component WARMABS fit, i.e. model 4. That is, we have included an additional emission component, 
which emits due to 'thermal' (i.e. collisional and recombination) processes rather than resonance 
scattering. This component is calculated using the photemis model, and is the emission analog 
of WARMABS. It is assumed to have the same ionization parameter, abundances, redshift, and 
turbulent velocity as the warm absorber component. We allow the normalization of the emission 



^http:/ /cxc. liarvard.edu/proposer/POG/html/chap8. html#fg:hetg-heglrf 
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component, which is proportional to its emission measure, to vary, along with the optical depth of 
the warm absorber. In doing so, we find the best fit to be that with zero thermal emission, and 
the statistical upper limit on the thermal emission component corresponds to an emission measure 
EM < 1.2 X 10^^ cm~^. We will discuss the implications of this in the next section. 

We have also examined the possibility that the absorption is produced in a plasma which is in 
coronal equilibrium instead of photoionized. This is what might be expected if the gas is heated 
mechanically, and if the outflow is due to thermal expansion of such a mechanically heated wind. 
This is done using the xspec analytic model hotabs, which calculates the absorption spectrum of 
partially ionized gas if the ionization is due to electron impact. The free parameter describing the 
ionization balance in this gas is the electron temperature. A key difference between a photoionized 
gas and a gas in coronal ionization equilibrium is that the ionization abundance distributions from 
photoionized gas have more overlap in parameter space than does a coronal gas. That is, at a given 
temperature each element in a coronal gas is most likely to exist in a pure ionization state, while at 
a given ionization parameter in a photoionized gas, each element is likely to have a mixture of two 
or more ionization states. In addition, the ions which can coexist in a coronal plasma all tend to 
have similar ionization potential. A consequence is that it is impossible for coronal equilibrium to 
allow the coexistence of H-like or He-like ions of elements with very different nuclear charge, eg. Fe 
XXV and S XVI. In contrast, a photoionized plasma at log(^)=4 does allow this. For this reason, 
coronal equilibrium models do not fit the GRO J 1655-40 spectra as well as photoionized models. 



3.5. Abundances 



We have also varied the elemental abundances and explored the limits for these allowed by 
the statistical criterion. In doing this, we rely on the criterion of Cash (19791, where a 99% 



confidence interval is defined by the values of the parameter which fall within Ax < 10 of the best 



fit value. We display these results graphically in figure 21 for the elements O - Ni. The abundances 



here are taken relative to the solar values of Grevesse et al. (19961 for abundant elements, andiAUen 



(19731 for the new elements added for this calculation. In all of our models we fix the abundance 



of Fe at 1 relative to solar. In Figure 21 the results of model 4 are shown in black, and the results 
of model 6 are shown in red. This shows that the models agree on the abundances of elements Sc 
through Co, such that the abundances of Cr, Mn and Co are all enhanced relative to Fe by at least 
50 %. For Ca model 4 predicts values greater than solar by 0.5-1 dex while model 6 is consistent 
with solar. As discussed above, we conclude that model 6 most nearly fits the overall properties of 
the spectrum, and in what follows we discuss the implications of the abundance pattern from this 
model. Limits on the abundance of O come from the O VIII Is - 5p line at 14.8 A and Is - 6p at 
14.64 A. This results in Yq values in the range 0.2 - 1.5. 



The abundance patterns shown in figure 21 are crudely consistent between models 4 and 6, 
and show enhanced abundances of Cr and Mn relative to Fe. Results for Ca-V are ambiguous, and 
suggest no strong enhancement, while the abundances of Na-Cl are smaller than the solar ratio. 
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relative to Fe. 



The overabundances of Fe-group elements strongly suggest that the observed matter has been 
subject to high-temperature burning conditions. For this reason, we have compared the observed 
abundances against nucleosynthesis models in massive stars, and more particularly resulting from 
the hydrostatic C-, O- and Si-burning phases. We use the presupernova nucleosynthesis yields 



calculated in model stars of 15, 20, and 25 solar masses with solar metallicity (Limongi et al 



20001. Each of the C-, O- or Si-burning zone has been weighted by the relative masses needed to 



best reproduce the observations. The relative mass ratio (C:0:Si) obtained is (6:1:1) for the three 



model stars. As shown in Fig. 22 the predicted abundances (relative to solar) are in rather close 
agreement with the observed pattern for models 4 and 6. 

As far as the light elements are concerned, the agreement between theory and observation is 
rather good. The O overabundance predicted reaches 1.3 (1.8 for the 20 Mq star) in agreement with 
the 1.5 upper limit determined in the 2-component model. Note, however, that the O abundance is 
not well constrained by the observed spectrum owing to strong interstellar absorption. Furthermore, 
we consider it possible that the low Z elements are affected by a separate ionizing continuum 
component. This could suppress absorption from O even if the abundance were greater than solar, 
although this is contrived given the absence of observations of such radiation. We also point out 



the possibility that the Ne and O lines are at least partly of interstellar origin (e.g. Juett et al. 
(2004|). This would require a coincidence in velocity between the GRO J1655 outflow and the 



intervening gas. It would also place more severe constraints on the nucleosynthetic models, since 
it would decrease the inferred abundances of these elements intrinsic to the source. 

Concerning the underproduction of Na, the disagreement may be due to uncertainties affecting 
the nuclear reaction rates or the thermodynamical conditions in the combustion zones. For the 
heavier elements above Ca, the agreement is still good though discrepancies can be seen in particular 
for Ti, V, and Co. It should be mentioned here that, in the nucleosynthesis simulations, all unstable 
nuclei produced have been assumed to /3-decay (except the long-lived ^^Mn). 



It is interesting to compare with the results of the optical abundance determination by |Israelian 
et al.| ( [l999 1. These authors find evidence for enhanced O/H, Mg/H, Si/H, S/H relative to solar, 
but Fe/H and Cr/H are approximately solar. This differs from our results for models 4 and 6, 



insofar as they can be compared, since Israelian et al. (19991 do not measure abundances for Mn 



Co, and Ni. The conclusions of Israelian et al. (1999) have not been confirmed in an investigation 



by Foellmi et al. (2007 1 . This also underscores the fact that the abundances we measure are relative 



to each other, and all are for elements heavier than O. We have assumed that iron is solar and quote 
other abundances relative to that, but we have no constraints on abundances for light elements, or 
for any abundances relative to hydrogen. 
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4. Discussion 

Other implications of our models include the fact that the curves of growth, in the absence of 
partial covering, indicate a velocity structure which has a small turbulent velocity, ~50 km s~^ at 
the same time as a larger bulk velocity, ~400 km s^^ . Since the radial velocity of the GROJ1655- 
40 system is 141^1 km s~^, this implies an outflow consisting of material which is cold or still 
compared with the bulk flow. This is unusual when compared with flows which are well studied in 
stars and non-thermal outflows from compact objects. 

All of our models fit to the iron K region by assuming the existence of a separate, higher 
ionization component which is primarily responsible for producing the Fe XXVI La line. This 
is because essentially all the other lines in the spectrum are consistent with a single ionization 
parameter and outflow velocity. Models 4 and 6 fit the Fe XXVI La line using an ad hoc Gaussian, 
while model 5 fits it using a separate warmabs component. Although the latter approach is more 
physically consistent, it does not fit the feature as well because for this model we force the two 
components to have the same outflow velocity. 

In either case, we find that the majority of lines fit to a single ionization component, log(^)~4. 



In addition to this, we have the observed luminosity L ~ 5 x 10'^^ erg s~^ (Miller et al. 20081 and 
density constraints from the Fe XXII metastable line implying n > lO^^cm"^. Taken together, 
these determine the location of the absorber: R = y^L/n/^ ~ lO^Lg^ n^^'^'^^^'''^ cm, and its size: 
AR = N/n ~ 10^A'^24'T'i5^ cm. In these equations L37 is the luminosity in units of 10^^ erg s~^, 
ni5 is the density in units of 10^^ cm~^ and .^4 = ^/lO^. That is, for plausible values, R <7 x 10^ 
cm, AR < 10^^ cm. This location can be compared with the 'Compton radius' within which 



photoionized gas cannot escape the gravity of the black hole (Begelman et al. 19831, which is 
Rjc = lO^'^(M/M0)T^g, where Tjcs is the temperature of the photoionized gas in units of lO^K, 
and for the GRO J1655-40 spectrum this value is <0.03 and the mass of the black hole is 5 - 8 Mq. 
If so, Rjc — 2 X 10^^ cm. Thus the inferred wind location is well within the Compton radius and 
is inconsistent with an outflow driven by thermal expansion, even though weak flows are possible 



at ~ O.lRjc as stated by Woods et al. (1996). This conclusion is consistent with that of Miller et 



al. (20081 



On the other hand, we note that the virial radius corresponding to the outflow speed we 
measure is 1.3 xlO^^ cm, considerably greater than the position we infer from the X-ray ionization 
balance. This, together with the fact that the outflow speed is comparable to the orbital speed in 
the system raise the possibility that the outflow could be associated with the secondary star or a 
region of the binary which is at comparable distance from the black hole. If so, then the radius we 
infer from the ionization balance arguments is an underestimate. This is difficult to understand in 
view of the density constraints from the Fe XXII lines; it would require that the metastable level 
of Fe XXII were populated by radiative pumping rather than by collisions. This, in turn, would 
require that the intensity of the radiation in the ~ lOOA region to be close to LTE, and there 
is no known source for such radiation far from the black hole. In addition, the observed outflow 
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speed does not appear to vary during the observation ( Miller et al.||2008 |, which spanned ~20% of 



an orbital period, which argues against an origin associated with the companion star or accretion 
stream. Thus, we consider this possibility to be unlikely. 

Another astrophysical system which appears to show a comparable contrast between line tur- 
bulent width and outflow velocity are FU Orionis stars. In these systems the terminal velocity of 
the wind is 300-400 km/s, and the rotational broadening (seen in absorption) of the wind is about 
50 km/s ( [Calvet et al.|[T993 Hartmann Sz Calvet||1995 |. The intrinsic turbulent velocity needed to 



produce the lines is likely much smaller. If so, the contrast between outflow and turbulent veloci- 
ties is comparable to the contrast between the observed width and the virial speed at the inferred 
position for the GRO J1655-40 outflow. 

We note also that the flow timescale is tdyn = R/vout ~ 27s. This can be compared with 
the recombination timescale tree = (na)^^ ~ 10^'^s, showing that the assumption of ionization 
equilibrium is valid for these simple assumptions about the density and location of the outflow. 
This does not change qualitatively if the absorber is located at the virial radius. 

We can also explore the implications of the emission measure upper limit derived in the previous 
section, EM < 1.2 x 10^^ cm~'^. The emission measure expected from a constant density shell with 
size derived from the ionization parameter, density and thickness is EM = 10^'^O-L37A''24^^^cm~^, 
where O is the solid angle of the warm absorber /emitter. So we infer 17 < 0.12Lj^y A^^^^4 or 
Q < 0.024 for the most likely value of -L37. This can be compared with the constraint derived by 



Miller et al. (20061 which is $7 < 1.4 based on the observed limit on the Fe XXIV 2p - 3s emission 
line at 11.43 A. We can also then estimate the mass loss rate in the outflow: M= ^IB^nmuVout = 
2 X lQ^^N2^vj gm s""*^ where vj is the outflow speed in units of 100 km s~^; so M~ 8 x 10"^^ gm 
s^^ for our flts. This is small compared with the mass accretion rate required to fuel the X-ray 
Mace — L/{rjc^) ~ 6 X 10^'' gm The Miller et al. (20061 limit on the solid angle allows 



a considerably larger mass loss rate, <10% of the accretion rate. 

Another conclusion of our work is the existence of two velocity components, one with ionization 
parameter log(^)=4 and speed ?;=375 km s^^ relative to the Earth, and the other with higher speed 
and much greater ionization parameter. We flnd for this second component a speed ~1400 km s~^, 
and infer that the ionization parameter be high enough that Fe XXVI and Ni XXVIII are the only 
ions with appreciable abundance, say log(^) > 6. Then the mass flux is M= lO^^^g^^ug^ gm s~^, 
where we have assumed a source luminosity 5 x 10^^ erg s~^ and fg is the speed in units of 1000 km 
s~^. It is more difficult to constrain the fractional solid angle of this component than it is for the 
lower ionization component because any emission may be masked by the lower velocity component, 
and we have not attempted to test this. It is clear that this component can carry more mass than 
the low ionization component if its solid angle exceeds 0.02. However, the radial location of this 
component cannot be reliably established, since the features in the spectrum which provide density 
constraints do not share the outflow velocity of this component. 



The results found here agree qualitatively with those of Miller et al. (20081 in the sense that 
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the ionization parameter and density of the flow containing most of the lines imply a wind location 
which is too close to the black hole to be easily explained by a driving thermal mechanism. We 
favor a scenario in which there are two components to the flow, and the lower ionization component 
is somewhat less ionized and the higher ionization component is more ionized than the single 
component found by Miller et al. (2008). Our spectral fitting allows us to constrain elemental 



abundances, and we find enhanced abundances for several elements in the iron peak. We do not 
find a pattern of systematic enhancements for a-capture elements. 

This work was funded in part by a grant from the Chandra theory program. We thank the 
referee. Frits Paerels, for many constructive comments. 



A. Appendix: Photoionization Models 

In order to model the spectrum of GRO J1655-40 we have made modifications and enhance- 
ments of the XSTAR 1^ photoionization code ( jKallman and Bautista|[2001 1. This code calculates the 
transfer of X-rays and other ionizing radiation, and the ionization balance, opacity and reprocessed 
emission from gas under a variety of physical conditions. It contains a relatively complete and 
up-to-date collection of atomic cross sections, rate coefficients and atomic energy levels. The code 
and atomic database, along with the 'wARMABS/PHOTEMis' analytic models for xspec, are freely 
available, distributed as part of the ftools package, and have been widely used in interpreting X-ray 
spectra. 

The standard XSTAR distribution includes all abundant even-Z elements (plus N), i.e.: H, He, 
C, N, O, Ne, Mg, Si, S, Ar, Ca, Fe and Ni. In order to model GRO J1655-40 we have added ah the 
other elements with Z < 30, with the exception of Li, Be, and B. The atomic data needed to do 
this, in addition to the resulting ionization balance for an optically thin gas, are discussed in this 
appendix. These data include: energy levels, line wavelengths, oscillator strengths, photoionization 
cross sections, recombination rate coefficients, and electron impact excitation and ionization rates. 
Distorted wave dielectronic and radiative recombination rate coefficients, photoionization cross 



sections, and photoexcitation-autoionization rates have been calculated (Badnell et al. 2003) for 



ions of these elements belonging to the H-like through Na-like isoelectronic sequences and are 



also available in the ORNL ADAS database (Schultz, et al. 2000). Concerning the other atomic 



parameters, a survey of atomic databases (Kallman and Palmeri 2007) reveals that they are far 



from complete for any of these elements. Therefore we rely heavily on hydrogenic scaling for many 
of these quantities. 

In order to model the spectrum we have modified the xspec analytic model warmabs in order 
to account for the new elements. WARMABS uses the optically thin ionization balance calculated 



^http:/ /heasarc. gsfc.nasa.gov/docs/software/xstar/xstar. html 
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by XSTAR, along with the same atomic database as xstar, to calculate the synthetic spectrum 
expected for a gas seen in transmission with a given ionization parameter, column density, turbulent 
velocity and abundance set. This takes into account all of the absorption processes, i.e. lines and 
photoionization cross sections, included in xstar. A complementary model, photemis, calculates 
the emission spectrum expected from such a gas, but that is not used in the calculations presented 
here. 

Rates and atomic structure for elements H, He, C, N, O, Ne, Mg, Si, S, Ar, Ca, Fe and Ni are the 
same as those described in ( Kallman and Bautista|2001 Bautista and Kallman|2001 1, and available 
in version 2.11 of the XSTAR code ( http://heasarc.gsfc.nasa.gov/docs/software/xstar/xstar.htmll ). 
For the odd-Z elements, and the iron peak elements not included in the release version of XSTAR, 
we use atomic structure scaled according to the following prescription: for the H- and He-like 
isoelectronic sequences we scale all energies and rates according to that expected for hydrogenic 
quantities: energies scale according to Z^, photoionization cross sections scale according to Z~'^, 
transition probabilities scale according to Z^. For these isosequences we scale the full atomic 
structure and rate set from xstar, including all levels through n = 6 for hydrogenic and through 
n = 5 for helium- like ions. For other isoelectronic sequences we adopt hydrogenic scaling of a highly 
simplified hydrogen-like ion, consisting of two spectroscopic levels. Is and 2p, along with a super- 
level and continuum. The energies, transition probabilities, and other rates are scaled according 
to a hydrogenic prescription with an effective Z = ^ E'j/j/lS.GeV, where is the first ionization 



potential taken from Allen (1973). Following this, we correct the energies of the n=2 levels using 



the evaluated wavelengths of the 2-1 transitions from the NIST database. These are available for 
F VIII, Al XII, Al XIII, Sc XX, Sc XXI, Ti XXI, Ti XXII, V XXII, V XXIII, Cr XXIII, Cr XXIV, 
Mn XXIV, Mn XXV, Co XXVI, Co XXVII, Cu XXVIII, Cu XXIX. We do not have evaluated 
wavelengths for the K lines of Na, P, CI, and K. 

For ions not previously included in xstar, ground state photoionization cross sections are 



taken from the calculations of Verner &; Yakovlev (1995). 



Recombination rate coefficients were taken from the calculations reported in the series of the 
papers by Badnell and coworkers ( Badnell et al.|[2003 | and available from the website]^ These are 
available for all elements He - Zn and for isoelectronic sequences H ~ Mg-like for both dielectronic 
(DR) and radiative (RR) recombination. For isoelectronic sequences Al-like through Ni-like, we 
make use of rates from [Aldrovandi and Pequignot (19731; Arnaud & Raymond (1992); Shull and 



Van Steenberg (1982), which cover ions of Si, S, Ar, Ca, Fe, and Ni in these isosequences. For ions 



of other elements in these isosequences we interpolate along isoelectronic sequence, xstar does 
not directly use the total recombination rate, but rather calculates rates onto a set of spectroscopic 
levels, typically with principle quantum numbers n <6, using photoionization cross sections and 
the Milne relation. Then it calculates a photoionization cross section for one or more fictitious 
superlevels such that the total recombination rate for the superlevel(s) plus the spectroscopic levels 



''http:/ /amdpp. phys.strath.ac.uk/tamoc/DATA/ 
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adds to the desired total rate taken from one of the above compilations. The super levels are 
generally assumed to decay directly to ground without the emission of any observable cascade 
radiation. The exception is the decay of the H- and He- isoelectronic sequences, for which we have 
explicitly calculated the decay of the superlevels to the spectroscopic levels using a full cascade 
calculation (Kallman and Bautista 2001; Bautista and KallmanpOOl |. 



A.l. Ionization Balance 

A key input affecting the ionization balance is the spectral energy distribution (SED) assumed 



for the source. We take the spectrum adopted by Miller et al. (20081, consisting of a disk black 
body with inner temperature kT=1.35 keV, plus a power law with photon index T = 3.54. The 
normalizations are chosen such that the low-energy cross-over between the two components is at 
1 keV. Below this energy we assume that the spectrum is flat (in photons). We assume that the 
spectrum incident on the gas in the GROJ1655 system is not affected by interstellar absorption. 
This spectrum is much steeper in the X-ray band than the conventional T ~2 power law applied to 
AGN. The ionization balance is correspondingly quite different in the sense that the mean ionization 
is lower at a given value of the ionization parameter in the steep spectrum case. 

Our choice of spectrum is deficient in photons between ~1 and 100 Ry, relative to higher ener- 



gies. This leads the gas to be thermally unstable, as originally described by Buff Sz McCray (1974 1 . 
This manifests itself as a region of ionization parameter where the temperature and ionization 
balance of the gas can be multi-valued. In the case of the spectrum we have chosen, this occurs at 
an ionization parameter log(^) ~ 2. This is less than the range of ionization parameters which can 
produce the ions observed from GRO J 1655, and so it will not affect our results. In what follows 
we will not discuss this further. 



Figures 23]|26 show the ionization balance for our choice of incident spectrum, as a function of 
ionization parameter ^ = L/{nE?). In our case, we adopt cgs units, so has units erg cm s^^. It 
is clear from this figure that the dominant observed ions require ionization parameters log(^) >4, 
but that production of Li-like and B-like iron and Ni requires ionization parameter values at the 
low end of this range. 
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Fig. 16. — Element column densities obtained assuming linear curve of growth, cosmic element 
abundances, and unit ion fractions for all lines in Table [2j Each data point represents one line. The 
horizontal axis is the quantity Zjon where we distinguish between contributions of various ions by 
assigning Zjon = -^eiement + 1 — 0.1 * (ionstage), where {ionstage)=l for hydrogenic, 2 for He-like, 
etc. Element color coding is the same as in figure [T5| 
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Fig. 17. — Curve of growth for H-like and Li-like ions. Element color coding is the same as in 
figure 15 Dashed diagonal lines show simple proportionality behavior expected for unsaturated 
lines. Solid lines are linear regression fits to data with errors less than 0.1. 
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Fig. 18. — Ratio of column densities derived from linear curve of grow analysis for H-line to He-line 
ions, versus the nuclear charge. Element color coding is the same as in figure 15 
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Fig. 19. — Curve of growth for various choices of the velocity characterizing the Doppler broadening 
Vfurbi and for a damping parameter corresponding to that for the Si XIV La line, at a Doppler 
broadening velocity of 20 km s~^. 
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Fig. 20. — detail of spectrum in the 1.5 - 2 A region. Crosses are data, blue is model 4, red is 
model 5, green is the Fe XXVI La Gaussian contribution to model 4. 



- 51 - 




-1.0 



CD CO a)_ ._ CO o _ 1- ^ O 

OLLZZS<a3Q_C/DO<^OC/3H>OSLLO 



element 



Fig. 21. — Element abundances: model 4(black) and model 6 (red). Iog(abundance)=0 corresponds 
to solar Grevesse et al. (1996) values. 
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Fig. 22. — Comparison of observed and predicted overabundances. The symbols with errors bars 
correspond to the 1-component (black) and 2-component (red) models. The lines correspond to 
model calculations in a 15 Mq (solid), 20 M0 (dashed) and 25 Mq (dotted) star, as described in 
the text. 



- 53 - 




log© 



Fig. 23. — ionization balance: H-Ne 
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g. 24. — ionization balance: Na-Cl 
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Fig. 25. — ionization balance: Ar-Cr 
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